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Abstract 

Combining feebly interacting massive particle (FIMP) dark matter (DM) with scale invariance 
(SI) leads to extremely light FIMP (thus the FImP) with FImP miracle, i.e., the mass and relic 
generations of FImP DM share the same dynamics. In this paper we show that due to the lightness 
of FImP, it, especially for a scalar FImP, can easily accommodate large DM self-interaction. For a 
fermionic FImP, such as the sterile neutrino, self-interaction additionally requires a mediator which 
is another FImP, a scalar boson with mass either much lighter or heavier than the FImP DM. DM 
self-interaction opens a new window to observe FImP (miracle), which does not leave traces in the 
conventional DM searches. As an example, FImP can account for the offsets between the centroid 
of DM halo and stars of galaxies recently observed in the galaxy cluster Abel 3827. 
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The conventional weakly interacting massive particle (WIMP) paradigm for dark matter 
(DM) is being challenged by quite a few null results of direct and indirect DM searches. 
They are already probing the typical WIMP DM, and yield particularly strong bounds in 
the lighter DM region. On the other hand, viewing from DM relic density along, WIMP does 
not take advantage over FIMP, i.e., the feebly interacting massive particle [1, 2], Instead of 
the freeze-out dynamics of WIMP that is used to keep thermal equilibrium with the plasma, 
FIMP gains correct relic density Vth 2 ~ 0.1 through the freeze-in dynamics. It never enters 
the plasma but still arrives Vth 2 ~ 0.1 via slow thermal productions like thermal particles 
decay. An obvious merit of FIMP is that it, just as expected, leaves null results in the 
conventional DM detectors devised for WIMP DM. Moreover, as a competitor to the WIMP 
miracle, a miracle of FIMP arises from combining FIMP with scale invariance (SI) [3]. This 
classical symmetry may provide a way to address the hierarchy problem [4]. 

FIMP with SI is a nontrivial combination, which gives rise to several important conse¬ 
quences. Immediately, quite generically FIMP must be extremely light (thus dubbed FImP). 
The point is simple. In most of the SI schemes for generating the electroweak (EW) scale, 
spontaneously breaking of SI happens at the electroweak (EW) scale [5] or TeV scale [6], 
by means of a scalar field collectively denoted as <3? with vacuum expectation value (VEV) 
u = ($) < TeV. By virtue of SI, all particles including DM X should gain masses via 
coupling to scalar fields with non-vanishing VEV, for instance, to <f>. Schematically, we can 
write down terms for mass generation 

(X* 2 * 2 , i W V 2 4>, (1) 

with X assumed to be a real scalar and Majorana fermion, respectively. Because X is a 
FIMP, one has A 0 , 2 /^, 1 thus very light FIMP. The ensuing important consequence is the 

aforementioned FImP miracle, which now becomes obvious: (due to SI) the mass and relic 
generations of FIMP share the common dynamics. The final consequence is that dark parity 
by hand may be not necessary for FImP. It can be sufficiently long-lived even without an 
exact protective symmetry, because its decay width is greatly suppressed by light mass and 
moreover feeble couplings. 

FIMP barely produces observable signatures except for a decaying one [7] (but may leave 
hint in tensor-to-scalar ratio [8]). However, FImP does. Firstly, Since FImP, such as a sterile 
neutrino or Majoron [9], does not require a parity, it is well expected that it can decay into 
X —ray photon(s) l . Secondly, also the core of this paper, FImP is likely to have appreciable 
self-interaction 2 partially by virtue of its lightness; hence, that kind of FImP can be probed 

1 Neutrino and photon are almost massless particles in the SM, so probably they are the only kinematically 
accessible final states for FImP decay. Neutrino is hard to observe while X—ray line is a good observable. 

2 DM self-interaction was originally motivated to address the small scale problems [10]. 
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from its astrophysical effects, e.g., leading to a separation between the DM halo and stars 
of a galaxy which is moving through a region with large DM number density. Interestingly, 
such a phenomenon was reported recently [11], It was discovered in the galaxy cluster Abell 
3827, within which (in the inner 10 kpc core) four elliptical galaxies were observed and their 
DM halos were reconstructed, with at least one spatially offset from its stars by a distance 
of A = 1.62±g;tg kpc [11]. Such kind of offset can be explained by DM with self-interaction 
that leads to DM self-scattering rate per mass odm/^dm ~ (1.0 — 1.5)cm 2 /g [12], or in the 
particle physics unit: 

cr dm/wdm ~ (4.7 — 7.0) x 10 3 GeV 3 . (2) 


Noticeably, despite of some tension with the bullet galaxy cluster upper bound, this value 
is also indicated to solve the small scale problems [13]. Therefore, it is of great interest to 
explore FlrnP with self-interactions (see other attempts [14-16]). 

In the absence of velocity dependence, Eq. (2) is saying that the mass scale involved 
in DM scattering, mass of a mediator or DM itself, should be far below the weak scale. 
Actually, for a typical WIMP, a^u/m is expected to scale as 


ctdm 


1 A 2 


~ 10" 10 x I — 


A 

oTi 


/ lOOGeVY 

\ m J 


GeV 


-3 


( 3 ) 


m 32n m 3 

which is about 14 orders of magnitude smaller than the indicated value given by Eq. (2). 
By contrast, if the dark sector mass scale m ~ (P(MeV), we will get the correct order easily. 
But it immediately raises two questions: Is there any theoretical motivation for that light 
DM scale? And for that light DM how does it get correct relic density? Bare in mind that 
in the SM maybe only the photon and neutrino can be the final states of MeV scale DM 
annihilating, thus the second question concerns us. 

For our FImP with SI, the two questions are simultaneously addressed in a coherent way. 
In the light of FImP miracle, interactions from Eq. (1) are supposed to freeze-in dark matter. 
Concretely, it is the two-body decay $ —> XX that dominates the freeze-in processes. And 
then the final yield can be formulated as [2, 7] 

45 g* r(7/2)r(5/2)M P1 


Yx(oo) « 


1.667T A g*\fgZ 16 m\ 


T($ -> XX), 


( 4 ) 


with g<$, and m<$> the internal degrees of freedom and mass of <P, respectively. The parameters 
9* ( 9* ) are the effective number of degrees of freedom contributing to the entropy (energy) 
density at T ~ m$. Within the SM <?* ~ g* ~ 106. For multi mother particles contributing 
to freeze-in X, there is a summation over $. Eventually, with Yx(oo) one can express the 
FImP relic density as 


n x h 2 = 2.82 x 10 


5 / mDM \ 
\MeV / 


Yx(oo) 


( 5 ) 
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In the ideal FImP miracle case, it is proportional to X^ 2 and y^, for a real scalar and 
fermionic FImP respectively. For the typical scale of u and m$, which are not far from the 
TeV scale, Q x h 2 — 0.1 uniquely determines the feeble couplings A^ In the following 

two subsections, we will detail the scalar and fermionic FImP with large self-interactions. 


A. The scalar FImP 


A scalar FImP X = S can be easily realized in the scale invariant SM where only the 
Higgs doublet H radiatively obtains VEV 3 . And the resulting model embodies the ideal 
FImP miracle. The relevant part of the model is very simple (see relevant studies of real 
scalar as a FIMP without SI [18]): 

^S 2 \H\ 2 + jS 4 . ( 6 ) 

Note that here S is automatically stable as a consequence of SI [17]. This model leads 
to u = v = 246 GeV. After EW spontaneously breaking (EWSB), the first term solely 
determines DM mass and relic density. The FImP mass is ms = \JX s h/2v. Freezing-in of 
S is dominantly proceeding via the decay of the SM Higgs boson h into a pair of S, with 
decay width 

r (h -a ss) = —(7) 

y J 32n m h K J 


Then, in terms of Eq. (4) the relic density is estimated to be 


Vt x h 2 ~ 0.12 x 


A 


sh 


10 


- 10.5 


5/2 


v/mh 

2.0 


1(F 




( 8 ) 


Equating Q x h 2 with 0.11 one can fix the unique free parameter A s h and hence the mass of 
FImP, ms = 1.0 MeV. It was already obtained in the Appendix of our earlier work [3]. 

In most cases, the coupling constant A plays no roles in DM phenomenologies. However, 
this ignored parameter is potential to generate a large DM self-interaction, especially for the 
FImP scenario under consideration. The resulting FImP self-scattering rate in per unit DM 
mass is given by 


cdm 

m s 


1 A 2 
1287T m| 


~ 7.9 x 10 3 


/ MeV A 3 
\ rn s ) 



( 9 ) 


3 It is well known that this model fails in triggering successful clectroweak spontaneously breaking and 
then modifications are indispensable. But this is not of our concern here [5, 6, 19]. Our discussion is 
particularly suited for the modification where additional bosonic states are introduced to overcome the 
top quark. 
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The scattering is from the contact interaction of S, so the scattering rate involves only the 
DM scale. Given a light DM scale, the self-scattering rate easily becomes large as long as A 
is not very small. This general advantage for light scalar dark matters, which always allow 
a quartic self-coupling to generate significant self-interactions, has already been utilized by 
the early references [20]. We would like to stress again, not only the light DM scale but also 
correct DM relic density, which is fairly problematic for the MeV scale DM, are naturally 
and coherently achieved here by virtue of the FlmP miracle. 


B. The self-interacting femionic FlmP 

Although this example does not give an ideal FlmP miracle, it takes a theoretical ad¬ 
vantage, i.e., a Majorana FlmP DM candidate X = N is naturally predicted rather than 
introduced in the very low scale seesaw mechanism [21]. Its scale invariant version shows 
several theoretical merits [3]. Scalar singlets Si with non-vanishing VEVs are necessary in¬ 
gredients of the model, to generate Majorana mass for N. On the other hand, these singlets 
are also badly needed to implement hidden SI spontaneously breaking. At the same time, 
they alleviate the serious relic density problem of sterile neutrino DM through the freezing- 
in mechanism, admitting the FlmP miracle (not ideal, see reasons later). Both are not far 
from the TeV scale. 

Before heading toward the explanation to the not ideal miracle, here we brief the La- 
grangian and report some relevant conclusions. We refer to Ref. [3] for more details. Without 
imposing any symmetry by hand, the relevant Lagrangian takes a form of 

-C N = V{Si , H) + yjHN + jSiN 2 , ( 10 ) 

V(Si,H) is a generic scalar potential for the singlets and Higgs doublet, and its concrete 
form containing two singlets can be found in Ref. [3]. In practice, at least two singlets (here 
we use the minimal number) are needed to trigger EWSB and moreover accommodate a 
quite SM-like Higgs boson near 125 GeV [22], For later use, we denote the extra singlet-like 
Higgs bosons as H 2 and V, with the latter (also the lighter) one being the pseudo Goldstone 
boson of SI spontaneously breaking. 

The FlmP miracle is not as ideal as that of the scalar FlmP, mainly owing to the presence 
of multi singlets with VEVs. They provide multi sources for generations of FlmP mass and 
relic, whose numerical correlation hence is weakened. Even then, the miracle still holds in 
the sense of order of magnitude. Following Ref. [3], decays H a — y NN freeze-in DM and its 
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relic density can be parameterized as 


^DM ^ 2 


0.11 x 



mo m 
0.1 MeV 


3 ^lOTeVy 


/1000 GeV\ 

V m H a J 


10 3 \ 


(ii) 


where tudm = = JY A i(Si ). Compared to Eq. (8), the extra parameters fn a (also the 

undetermined masses mjjJ manifest the deviation from the ideal FImP miracle. They are 
model dependent, on the patters of singlets VEVs and as well their coupling to N. But in 
most cases they are order one numbers. Note that the fermionic FImP is favored to be in 
the much lighter region, around the sub-MeV scale or even below given that the singlets 
VEVs are not far above the TeV scale. The reason is nothing but that the fermion mass is 
proportional to the coupling constant instead of its square root like the scalar FImP. To 
maintain the coldness of the FImP DM, we take a conservative value Mn = 0.1 MeV in this 
paper. Probably, it can be lowered down substantially, on account of a mildly colder DM 
spectrum from freeze-in [3, 23] 4 . 

Now we turn our attention to self-interaction of N. Unlike the scalar FImP, large self¬ 
interactions are not a built-in part of the fermionic FImP. It calls for a light mediator, either 
a vector or scalar boson [24], Viewing from our Lagrangian Eq. (10), a light scalar boson, 
denoted as So, is a natural choice. So is also a FImP, but it has a sizable coupling to N via 
the Yukawa coupling C So D —XS 0 NN (in four components). With it, N can scatter with 
each other through s/t/u— channel exchanging So- To get a sufficiently large scattering rate, 
two options are of interest here. One is a heavy So with mass ms 0 much larger than Mn 
and the other one is the opposite. In what follows we discuss them case by case. 


1. Dark force 

If Sq is very light, it becomes a dark force mediator and DM self-scattering, in the non- 
relativistic limit, is described by the following attractive Yukawa potential 

V = e - ms ° r . (12) 

r 

with a\ = A 2 /47t. In different parameter space spanned by {otd, M]y,ms 0 ), the potential 
may induce different velocity-dependent DM self-scattering, and we refer to Ref. [13] for a 
comprehensive discussion. Here we focus on the simplest case, i.e., a\M N /ms 0 <C 1 such 
that the Born approximation holds. Then, the perturbative computation in ot\ from V leads 

4 Despite beyond the scope of this paper, it is of interest to investigate the cosmological implications of 
warm DM with self-interaction. As shown here, a fermionic FImP with miracle tends to be a warm DM. 
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to [13] 


^Born 
Cj t 


&'Ka\ 

M> 4 


*°g (i + C) 



(13) 


with £ v = M^v/mso- Note that in the small /ms 0 region such that £ c 1 one actually 
gets the velocity-independent approximation cr^ orn ~ 47 t a\M‘f 4 /m A So . If the self-interaction 
could leave observable effect at the cluster scale without spoiling the small scale structures, 
we should work in this limit 5 . As an estimation, we typically need parameters as 


a, 


Born 


M, 


= 7.9 x 10 3 


N 


M, 


N 


0.1 MeV 


0.002 MeV 


™s 0 J 




GeV 


-3 


(14) 


The resulted = 50v, which is indeed a small number for the typical velocity v < 10~ 3 . 

We make a comment on the fate of So. It is at the keV scale, assumed to gain mass as 
the scalar FImP in Eq. (6). At leading order, it can decay into a pair of neutrinos, induced 
by the tiny active-sterile neutrino mixing, < 10~ 10 . The resulting decay lifetime > 10 24 s is 
much longer than the cosmological timescale, so it survives as a relic today. But its energy 
fraction is negligible due to its lightness and small yield during freeze-in. 


2. Four-fermion interaction 


The other option is a heavy So, which actually leads to contact four-fermion interaction 
again, but here we make a direct calculation of the self-scattering without integrating out So- 
This process receives all s/t/u —channel contributions, and from them we get the following 
scattering rate per DM mass 


0dm 

M n 


3A 4 M n 
87r m 4 o 


~ 6.0 x 10 3 


( M N \ 

Vo.lMeVy 


/ MeV \ 4 

V m s 0 ) 



(15) 


So is not favored to be much heavier than the MeV scale due to two reasons. One is that A 
will become accordingly large, even larger than order 1,which is unpleasant at low energy. 
The other one is to prevent it from spoiling the FImP miracle. So here is like a frozen- 
in scalar considered in Ref. [25], because itself is a FIMP and moreover could produce N 
via decay So —> NN. Therefore, the heavier So means the larger yield of Sq thus larger 
contribution to N production. Let us estimate this from Eq. (8) which just parameterizes 
the relic density of FImP S with mass 1 MeV. As a naive estimation, one can get the N 
relic density inherited from So decay by multiplying Eq. (8) a factor /ms 0 ~ 0.1. Thus, 


5 Otherwise, self-scattering is over enhanced in the small scale system like dwarf galaxy with characteristic 
v ~ lOkm/s but is insufficient in the large scale system like cluster where v ~ 1000 km/s. 
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this contribution to the final relic density of N is subdominant and the FlrnP miracle is not 
significantly affected. However, as we take heavier and heavier So, it will quickly dominate 
over the direct freeze-in production of N 6 . 

To end this paper we make a summarize. FlrnP is a necessary result after the combination 
between FIMP and SI, which further creates a FlrnP miracle. We show that a large DM self¬ 
interaction can be easily accommodated for FlrnP due to its lightness. This is particularly 
true for a scalar FlrnP which always has a quartic self-coupling. While for a fermionic 
FlrnP one has to introduce a mediator which is another FlrnP, a scalar boson with mass 
either much lighter or heavier than the FlrnP DM. DM self-scattering opens a new window 
to observe FlrnP (miracle), which does not leave traces in the conventional DM searches. 
For instance, they are potential to explain the recently observed DM self-interaction in the 
galaxy cluster Abel 3827. 

In the late stage of this paper, we found that Ref. [1] basically already studied the model 
of SM extended by a real scalar FlrnP, aiming at solving small scale problems using the 
self-interaction of FlrnP; Moreover, it pointed out that the MeV mass scale is consistent a 
model with zero bare Higgs mass, which is nothing but the classical scale invariance in our 
FlrnP framework. 
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